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Abstract 

Transport coefficients can be obtained from 2-point correlators using the Kubo formulae. It lias 
been shown that the full leading order result for electrical conductivity and (QCD) shear viscosity 
is contained in the re-summed 2-point function that is obtained from the 3-loop 3PI re-summed 
effective action. The theory produces all leading order contributions without the necessity for power 
counting, and in this sense it provides a natural framework for the calculation. In this article we 
study the 4-loop 4PI effective action for a scalar theory with cubic and quartic interactions, with 
a non-vanishing field expectation value. We obtain a set of integral equations that determine the 
re-summed 2-point vertex function. A next-to-leading order contribution to the viscosity could be 
obtained from this set of coupled equations. 
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INTRODUCTION 



To date, little progress has been made on the calculation of transport coefficients beyond 
leading order. A next-to-leading order contribution to shear viscosity in pure 4 theory has 
been calculated in [1]. The authors take into account the shift in the soft scalar propagator 
pole due to the thermal mass. This produces a correction of order v^A relative to the leading 
order result. The momentum diffusion coefficient for a non-relativistic heavy quark has 
been calculated at next-to-leading order in [2, 3]. Hydrodynamic coefficients have been 
calculated from a kinetic theory approach using a second order gradient expansion, and 
working at leading order in the coupling [4]. In this paper we propose a new strategy to 
calculate transport coefficients beyond leading order using an n-particle irreducible (nPI) 
effective theory. 

An nPI effective theory is defined in terms of n functional arguments which correspond 
to a set of n-point functions that are determined self-consistently through a variational pro- 
cedure. This procedure re-sums certain classes of diagrams, and represents a re-organisation 
of perturbation theory (for recent reviews see [5, 6]). A well known example of a case in 
which selective re-summations play an important role in quantum field theory is the hard 
thermal loop theory, which includes screening effects that regulate infra-red divergences. 
nPI approximation schemes are of particular interest because they can be used to study far 
from equilibrium systems [7-12]. To date, numerical calculations have only been done for 
2PI theories where it has been shown that the convergence of perturbative approximations 
is improved (see [5, 6, 13] and references therein). 

Recently, it has been demonstrated that nPI effective theories provide a natural frame- 
work to organise the calculation of transport coefficients. A large Nf leading log calculation 
of QED conductivity and shear viscosity was done in [14] using 2PI effective theory. The 
electrical conductivity and QCD shear viscosity have been obtained at leading order from 
the 3-loop 3PI effective action [15-17]. The formalism produces integral equations for the 
vertex functions that re-sum the pinching and collinear singularities and produce the full 
leading order result. In this sense, the 3-loop 3PI effective action respects the natural or- 
ganisation of the calculation: it produces the full leading order result without the need for 
any kind of power counting arguments. Since power counting is notoriously difficult, and 
becomes increasingly complicated at higher orders, nPI effective theories could provide a 
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useful method to organise the calculation of transport coefficients beyond leading order. 

In spite of the complexity of the nPI effective action, a systematic expansion can be done 
in a self-consistent way. A self-consistently complete loop expansion of the effective action 
can be based on an equivalence hierarchy [18]: to 2-loop order, the infinite-PI effective action 
is equivalent to the 2PI effective action, to 3-loop order, the infinite-PI effective action is 
equivalent to the 3PI effective action, etc. We work with the 4-loop 4PI effective action 
for a scalar theory with cubic and quartic interactions, with non-vanishing field expectation 
value. The goal is to take a first step toward the calculation of next-to-leading order transport 
coefficients. We derive the integral equations that would produce the next-to-leading order 
contribution to the viscosity. We remark that the numerical solution of these equations is 
more difficult than anything that has been accomplished so far. 

In addition, there are unresolved issues that would arise if we attempted to extend the 
calculation to gauge theories [19, 20]. The renormalizability of a theory is related to the 
existence of symmetry constraints on the n-point functions. For nPI effective theories, 
symmetries and renormalizabilty are connected to the fact that proper n-point functions 
can be defined in more than one way. All definitions are completely equivalent for the 
exact theory, but they are not the same at finite approximation order. These issues are well 
understood for scalar theories and QED at the 2PI level. For scalar theories one can define a 
2-point function that satisfies Goldstone's theorem in the broken phase [21, 22]. For QED one 
can define n-point functions that obey traditional Ward identities [23, 24]. These symmetry 
constraints allow one to construct a complete renormalization that preserves the symmetries 
of the original theory [22, 25-28]. For nonabelian theories, the situation is more involved. 
It has been shown that at any order in the approximation scheme, the gauge dependence 
of the effective action always appears at higher approximation order [29, 30]. However, the 
gauge symmetries of the n-point functions are more complicated than for abelian theories, 
and renormalizability remains an open question. 

The definitions of the various n-point functions that we will use are given below. The 
symmetry properties of these functions, which are described in the corresponding paragraphs, 
have not been established for nonabelian gauge theories. 

(1) 1PI vertex functions are obtained by functional differentiation of the 1PI effective 
action. These functions satisfy the standard symmetry constraints, which reflect the sym- 
metries of the Lagrangian. For example: the 2-point function satisfies Goldstone's theorem 
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for scalar theories in the broken phase, and the Ward identity for QED (transversality in 
momentum space). 

(2) Variational n-point vertex functions are the functional arguments of the nPI effective 
action. Extremising with respect to these functions produces integral equations for the 
vertices called equations of motion. The vertices can be determined self-consistently by 
solving the equations of motion using a variational procedure. They do not satisfy standard 
symmetry constraints. 

(3) Re-summed n-point vertex functions are obtained by taking functional derivatives of 
the re-summed effective action, which results from substituting the self-consistent solutions 
for all variational n-point vertex functions with n > 2 into the nPI effective action. Re- 
summed n-point vertex functions have a geometrical interpretation in field space, and satisfy 
standard symmetry constraints. 

(4) Mixed n-point functions can be defined by differentiating self-consistent vertices, 
which define the extremum of the nPI effective action, with respect to the field expectation 
value. Integral equations for these vertex functions are obtained by functionally differentiat- 
ing the equations of motion of the nPI effective action with respect to the field expectation 
value (see, for example, Eq. (21)). For QED, it has been shown that these functions also 
satisfy standard symmetry constraints [23] . 

The complete set of integral equations which give the re-summed 2-point vertex function is 
derived in this paper. Using the Kubo formula that relates the shear viscosity to the 2-point 
correlator, a next-to-leading order contribution to the shear viscosity could be obtained from 
this set of coupled equations. The re-summed 2-point function is obtained by functionally 
differentiating the re-summed effective action. The resulting expression contains mixed 
vertex functions. These mixed vertex functions satisfy integral equations which contain 
variational vertex functions. The variational vertex functions are determined by additional 
integral equations. 

The paper is organised as follows. In Section we define our notation and give our result 
for the 4- loop 4PI effective action (some details of the calculation are given in Appendix ). 
In Section we calculate the re-summed 2-point vertex function. In Section we calculate 
the integral equations satisfied by the mixed vertex functions that appear in the re-summed 
2-point vertex function. In Section we calculate the integral equations satisfied by the 
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variational vertex functions that appear in the mixed vertex functions (that appear in the 
re-summed 2-point vertex function). 



THE 4PI EFFECTIVE ACTION 

4PI effective actions were introduced in Ref. [31, 32]. They were first discussed in the con- 
text of relativistic field theories in Ref. [33]. The 3-loop 4PI effective action was calculated 
in Refs. [18, 31, 32, 34]. The fourth Legendre transform has been studied in [35]. In this 
article we study a scalar theory with cubic and quartic interactions, with a non-vanishing 
field expectation value. We use the method of successive Legendre transformations [18, 31] 
and obtain the 4PI effective action to 4-loop order. The result is given in this section, and 
some details of the calculation are presented in Appendix . We use the following notation: 
D°, U°° and V° are the bare propagator and vertices, D is the variational propagator, U is 
the variational 3-point vertex, and V is the variational 4-point vertex. We use a compact- 
ified notation in which a single numerical subscript represents all space-time co-ordinates. 
For example: the field expectation value is written 0i := <f)(x), the variational propagator 
is written _D 12 := D(xi,x 2 ), the bare 4-point vertex is written V® 2M := V°(xi, x 3 , x 2 , x 4 ), 
etc. We also use an Einstein convention in which a repeated index implies an integration 
over space-time variables. We will define propagators and vertices with factors of % so that 
figures look as simple as possible: lines, and intersections of lines, correspond directly to 
propagators and vertices, with no additional factors of plus or minus i. Using this notation 
we write the classical action: 

SM = \<f>l [i (D° 12 )- 1 ] 2 - ^U&fafafa - ^123401020304 • (1) 

We define the effective classical vertex U® jk := U°° k + 4>iV® jkl and obtain: 



o<p 2 o<pi 

123 503^02^01 503 

V ■ ggj _^1°23_ S\DU<P))-\ (2) 

1234 504503502^01 504 504503 
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The 4PI effective action can be written: 
r[4>,D,U,V}= SM + ^TrLnZ^ 1 

+ Z -Tt [(D^M)- 1 (D 21 - D° 21 (<j>))] + r°[0, D, U, V] + T mt [D, U, V] . (3) 

The propagator D and vertices U and V are to be determined self-consistently from the 
equations of motion. The terms T°[(f), D,U,V] and T mt [D, U, V] contain all contributions to 
the effective action which have two or more loops. The first piece r°[0, D, U, V] includes all 
terms that contain bare vertices. The calculation of T mt [D, U, V] is lengthy, but straightfor- 
ward. We show some of the steps in Appendix . The result is shown in Fig. 1 1 . We define 
$ = z(r° + r mt ) which allows us to represent $ as a series of diagrams with all factors of i 
absorbed into the definitions of the propagators and vertices. We have given each diagram 
in the figure a name, so that we can refer to them individually later. 

* = i oo+i 0-A0 + s©-ffiO 

EIGHT EGG EGG BBALL BBALL 




HAIR MERCEDES LOOPY TARGET EYEBALL TWISTED 

FIG. 1: 4- loop 4PI effective action. 

The equations of motion are obtained from the stationarity of the action. There are 
four equations which are obtained by functionally differentiating with respect to the four 
functional arguments of the effective action: 

X,e { *,D,U,V } . (4) 

The equations obtained by varying with respect to {D, U, V} can be solved simultaneously 
for the self-consistent solutions which are functions of the field expectation value: D[<p], 
U[4>], V [</)]. Substituting these self-consistent solutions we obtain the resummed action, 
which depends only on the expectation value of the field: 

f[0] = r[0, %>],%],%]]. (5) 



1 Figures in this paper are drawn using Jaxodraw. 
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In the future we will write T and T without their arguments. 



2-POINT VERTEX FUNCTIONS 

We start by defining some mixed vertex functions using the same notation as (2): 

O - vR - 5 ^ 123 (R\ 

"123 — Tl — ' ^1234 — — T" — • [0) 

These vertices are shown in Fig. 2. Legs that correspond to functional differentiation with 
respect to the expectation value of the field are called 'external.' These legs are distinguished 
by an arrow. For both Q and \1/ the last index is always assigned to the external leg. In 
order to simplify the form of the equations we will write Uo as Qq and Vq as \l/o throughout 
the next two sections. 



FIG. 2: The vertices Q and 

An additional useful relation can be obtained from the identity D^D^ = <5i2- Differentiat- 
ing with respect to <ft and using (6) gives: 

TT~ = -^H'-C > 22'^1'2'3 • (7) 

d(p 3 

The re-summed propagator is defined as: 

The re-summed 2-point vertex function, or the re-summed self energy, is extracted from the 
re-summed propagator using: 

( J DS t )- 1 = ( J D? 2 (0))- 1 -n» t . (9) 

We can derive an expression for the re-summed 2-point function as a function of the vertices 
in (6) by taking derivatives of the modified effective action and using the chain rule. We use 



the notation X; to indicate one of the set of functional variables X := {D, U, V} and Xi to 
indicate one of the set of self-consistent solutions X := {-D[0], U[<fi], V[<j>]}. We obtain: 



5 2 r 



+ 



2^01 \X 

5 2 T 



5T 



5X, 
SXi 



5 2 X 



X < 



(10) 



SXiSipi \x 50 2 



s 2 r 



5Xj5X,: 



5Xi SXj 



x i 



The last term in the first line is identically zero (see Eq. (4)). The expression can be further 
simplified by using the set of equations obtained by differentiating the equations of motion: 



_6_ 

502 



r ST 




I SXi 


X. 







S 2 T 



SXiS(p2 



5 2 T 



5Xj5Xi 



X 502 



0. 



(11) 



Using this constraint (10) becomes: 

i(Dtf)- 1 = 
Expanding the sum we have: 
i(Dtf)- 1 



5 2 T 



5 2 T 



5 2 T 



50 2 50i 



+ 



5 2 T 



5D 



M 



+ 



X (502 

5 2 T 



X 502 (5^345(501 



5U: 



345 



X 502 



(12) 



(13) 



5D 34 50i 

The last term in the sum does not contribute because the derivative 5 2 r/(5W>0) is identically 
zero. Using (2), (3), (6) and (7) we obtain: 



5 2 r 



5D 43 50i 

5 2 r 



SD. 



5U 5 4 3 5(f)i 



X 002 
5^543 



— — — 2(^0)341 • (-D33'-D44'^3'4'2) , — f2 + 2 — ^ 



X 502 



"g^3'4'5'l ' (^3'3-D4'4-D5'5^ / 3452) • 



(14) 



In the first line of Eq. (14), the two contributions on the right side come from the classical 
and 1-loop part of T (the first three terms on the right side of (3)). In the second line we 
have used 5/5-D 43 (5 < l ) /50i) = 5/5D43(5$egg o /50i) and defined the vertex fl' : 



5 2 T 



1 5 2 $ 

341 5D 4 350! 



2(^o)341 



'5L> 43 50i 2" d4i SD^Sfa ' 2 V " U '" 1 ' (15) 
where the term comes from the 1-loop part of T. The vertex Q' is shown in the first 
two graphs on the right side of the first line of Fig. 4. In the third line of Eq. (14) we have 
used: 

5 2 r 



'5U, 



543"V1 



^ Qggo _ 1 vrrO n n n 

0^543(501 6 



(16) 



If the future we will use the shorthand notation: 



DwDwDs&Vmm = (DDD^) 3 452 ■ (17) 
We extract U\f from (9), (13) and (14): 

nS?= ^234^3 + I (n' ) ul (DDQ), 32 + ^° U51 (DDm) 5432 . (18) 

The result agrees with the Schwinger- Dyson equation for the 2-point vertex function [36, 37] 
and is shown in Fig. 3. 



FIG. 3: The re-summed 2-point vertex function. 

The variational 2-point function satisfies an integral equation obtained from the equation 
of motion (4): 

2-^ = -+ {D 12 )-^{DUct>))^-Yl 12 . (19) 

By rearranging terms one can show that the variational 2-point vertex function satisfies the 
same Schwinger-Dyson integral equation as the resummed 2-point vertex function (Eq. (18) 
and Fig. 3), with the vertices fl and \I/ replaced by U and V respectively. 

INTEGRAL EQUATIONS FOR MIXED VERTEX FUNCTIONS 

In this section we derive the integral equations for the mixed vertex functions Q and \I/ 
which appear in the re-summed 2-point vertex function n ext (see Fig. 3). In addition to the 
vertices defined in (6), we will need to define the 5-point vertex: 

e 12345 = • (20) 

As in Eq. (6), the last index is always assigned to the external leg. This vertex is shown 
in the left side of the last line of Fig. 4. The integral equations satisfied by the vertices fl, 



\I/ and are obtained by taking functional derivatives with respect to the field expectation 
value of the appropriate equations of motion (4): 



<5 ' 


5T 




(509 




X. 


J_ 


' 5T 




509 


_5Ui'2'3> 


JL 


' 5T 


509 


_<5Vi'2'3'4 



= 0, 

n _1 n _1 f)~ l 



= 



D -1 D -1 

1J V\ 1J 2'2 1J 3'3 1J 



0. 



(21) 
(22) 
(23) 



The subscript X indicates that all self-consistent solutions are substituted. The numerical 
factors 2!, 3! and 4! are inserted for later convenience, as explained below. The index '9' on 
the derivative is chosen so that this index remains the last (highest) after using the chain 
rule (see Eqs. (24), (25), (27)). In Eq. (22) we multiply by three inverse propagators to 
remove the external legs that appear when we take the derivative with respect to U, and in 
Eq. (23) we multiply by four inverse propagators to remove the external legs that appear 
when we take the derivative with respect to V. 

The next step is to expand the derivatives in (21), (22) and (23) using the chain rule. We 
start with (21). For the 1-loop part of the effective action, we take the derivatives explicitly. 
Using (6) we obtain (including the factor of 2 that is introduced in (21) for this purpose): 
—^129 + (^0)129- Rearranging we have: 



^129— (^0)129 + 2 



5D. 



31 



<5 2 $ 



5D^5Di2 



x. 



+2- 



5U: 



345 



c5 2 $ 



5U 3i5 5D 



12 



+ 2- 



5V 



3456 



(5 2 $ 



5^3456 5£ 



12 



(24) 



Now we consider (22). The 1-loop part of the effective action gives zero contribution, since 
it does not depend on the vertex U. Using (6) the result can be written: 



1239" 



+3! 
+3! 



o 

1239 

su 456 



+ 3! 



5D 



45 



5$ 



5 



9 (5-D 45 L 

5$' 



SU V 



2'3' 



D^ 1 D -1 D -1 

xy l'l xy 2'2- Ly 3'3 



£09 <5£7 4 56 

<5U 4567 5 



dUi'2'3' 

5$ 



X 



4567 



D -1 D -1 D -1 

1^2' 2^3' 3 

A-l A-l f)-l 
xy l'l xy 2'2- Ly 3'3 



(25) 



_8Ui'2>3' 

The notation <3>' indicates that the diagram labelled EGG in Fig. 1 has been removed. This 
diagram, when multiplied by the factor 3!, produces the term — x I / i239- In addition, we have 
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used: 



5 2 $ 



oUi'2'3' 



D^ 1 D^ 1 

1J V\ iJ 2'2 1J 3' r i 



X 



= 3!— 


~$ 2 $EGG 


d0g 


(5L r i'2'3 / 



D -1 



^?239 • (26) 



Now we look at (23). The 1-loop part of the effective action gives zero contribution, since 
it does not depend on the vertex V. Using (6) the result can be written: 



©12349— 4! 

+4! 
+4! 



SD m 5 



SD 56 
5U 567 5 



5$ 



5Vi'2'3'4' 
5$ 



D -1 D _1 D^ 1 D -1 



X 



5( p9 5U 567 
5V 5678 S 



5Vi'2'3'4' 
5$" 



D -1 D _1 D^ 1 

1J V\ 1J 2'2 1J 3'3 1J <Vk 



^9 <5V £ 



5678 



SVi'2'3'4' 



fj-1 7-)-l 7-)-l 



(27) 



The notation <3>" indicates that the diagram labelled BBALL in Fig. 1 has been removed. 
This diagram, when multiplied by the factor 4!, produces the term — 612349. 
To simplify the resulting expressions, we introduce the following notation: 

£ 2 $ 



1 1 



C34;12 , 



5D 34 SD 12 2! 2! 
5 5<$> 11 

T77 T7^— = ^7^7^33'-D44'-D55'C3'4'5';12 , 

dt/345 <)L>i2 2! 3! 

5 5$ 1 1 1 
777 — xri = "77r~T7-^ 33 ' Du'D 55 'D 6& C3>4> 5 r 6 i.i2 + — <5i3-D44'-D 5 5'-D 66 'iV4/ 5 ' 6 / ; 2 , 



5 



6D 12 

5 

5U i5e 
5 

5^4567 
5 



5$ 



5Uy4' 5 ' 

5$' 



jj-i jj-i jj-i 



d~ x n -1 d~ 1 

±y l'l ±y 2'2- Ly 3'3 



#56 

^567 
5 



SUi'2'3' 
_5Uv2'3' 

5$ 



t-)-1 7-)-l 



- 1 1 C 

= -2j-2j-D44'-D55'-D66'C , 4'5'6' ; i23 + ^-^u-Dss'-tW^W^s , 



Vi'2'3'4' 

Vi'2'3'4' 
5$ 



d -1 d -1 

- ty l'l Xy 2'2- Ly 3'3- Ly 4'4 



D -1 



5678 



V v 



2'3'4' 



- D -1 D -1 D -1 



- 1 1 C 

~ ^j^j^56;1234 , 

= -^-7^13-D66'i}77'jV6'7';234 , 
111 

— -^-^^15^26^77' D S8 iN 7 r S i. 3 4 . 



(28) 



For each vertex, the semi-colon divides legs that attach to the left and right sides of the 
diagram. Using this notation and Eqs. (6) and (20), the equations in (24), (25) and (27) 
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can be written: 



^129 = (^o)l29 + 2^349-033' -D44'C3'4';12 

+ -^^3459-^33' -D44'-D55'C3'4'5';12 _|_ Q^ m _ D 32,' D iA i D 551 D C y y y §1 -^2 , 

^1239 = ^1239 + 2^ 4 59^ ) 44'-D55'C4 / 5';123 



(29) 



4569 



1 3 

— D 4 4/D 5 5/Z) 66 'C , 4/5' 6 / ; i23 + —^liD^iD m iN^ & . 2 -i 



+Q45679-^j^l4-D55'-D66'-D77'^5'6'7';23 , 

1 4 
©12349= -^569-D55' -D66'C5'6';1234 + _ 27^5679^15-D66'-t ) 77'^V 6 / 7';234 



4! 1 

+ 2^27 2 ©56789^15^26-077' -^88' N V8 ,. 



34 ■ 



(30) 



(31) 



Eqs. (29), (30) and (31) are represented diagrammatically in Fig. 4. 








22! 2 




FIG. 4: The integral equations for the vertices Q, \1/ and O. The grey boxes are the vertices 
listed in Eq. (32). 



We note that in Eqs. (29), (30) and (31) and Fig. 4 we have combined terms that 
correspond to permutations of external legs. For example, the factor 3 multiplying the 
fourth term on the right hand side of (30) indicates that three terms have been combined 
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as indicated in Eq. (32): 

2^ / 4519-D44'-D55'^4'5';23 + ^ ^4529-^44' -^55' N# 5'; 13 

1 3 

+ 2^4539-D44'-D55'^4'5';12 => ^ 4<o3$D Ai i D^i N w -i 2 • 

We introduce a shorthand notation for the vertices defined in (28) by suppressing indices. 
We write: 

^ij-yxy ■— ^22 j ^ijk;xy • — ^32 > ^ij;xyz • — ( - / 23 > ^ijk\xyz • — ^33 > ^vjkl;xy •— ^42 j , > 

(32) 

Cij-xyzw '■— C*24 , ^ij,xy '■= -^22 , Nijj^xy '■= -^32 , ^ij.xyz '■= N23 ■ 

These vertices are calculated by substituting $ (as shown in Fig. 1) into (28). In each case, 
the symmetric partner of a vertex is obtained by inverting the legs. For example, C23 is the 
right /left inverse of C32. To illustrate the procedure we show below the contributions to 
each of the vertices from the TARGET diagram. 



We obtain the vertex C22 by differentiating twice with respect to D. This corresponds 
to removing two propagators. Different contributions to C22 arise depending on which two 
propagators are removed. The result is shown in Fig. 5. 



C22 [TARGET] = 




FIG. 5: Contributions to the vertex C22 from the TARGET graph. 



The first graph represents 4 different diagrams that have all been combined in the figure to 
simplify the notation. These 4 diagrams correspond to the 2! 2! = 4 permutations of the 
legs on the left and right sides of the figure. The second graph contains a factor 2 instead 
of 4 because an additional factor of 1/2 is contributed by the symmetry factor of the graph. 
This symmetry factor comes from the fact that the diagram is symmetric under interchange 
of the two internal vertices. For the third diagram the factor is 2! 2! • |, where the | comes 
because permuting the left legs is equivalent to permuting the right legs (one can see this 
immediately from the fact that the diagram is symmetric when inverted top-to-bottom). 
The factor for the fourth and fifth graphs is obtained in the same way as for the first graph. 
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Next we consider the vertex C32. First, we differentiate with respect to U and remove 
the attached legs. Then we differentiate with respect to D. The results are shown in Fig. 6. 



C32 [TARGET] = 




FIG. 6: Contributions to the vertex C32 from the TARGET graph. 



The second and third graphs carry a factor 3! 2! = 12 that corresponds to the permutations 
of legs on the left and right side of the diagram. The first graph has an additional factor of 
1/2 from the symmetry factor of the graph, which comes from the fact that the diagram is 
symmetric under interchange of the two internal vertices. 

Next we calculate the contributions to the vertices C33 and N 2 2- We differentiate with 
respect to U and remove the attached legs, and then differentiate with respect to U again. 
In general, these derivatives produce two classes of terms: for some terms there is a delta 
function which ties together two of the legs (as in the term on the right side of the last line 
in Eq. (28)), and for some terms this delta function does not appear. The results are shown 
in Fie;. 7. 



C33 [TARGET] = 18 



N 22 [TARGET] = 4 




FIG. 7: Contributions to the vertices C33 and N22 from the TARGET graph. 



The first graph carries a factor 3! 3! • \ = 18 which comes from the possible permutations 
of legs on the left and right side of the diagram, with the \ accounting for the fact that the 
diagram is symmetric when inverted top-to-bottom. The second graph has a factor 2! 2! = 
4 from the permutations of the left and right legs. 

Next we calculate the contributions to the vertex C42. We differentiate with respect to 
D, and then differentiate with respect to V. The result is shown in Fig. 8. 
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C 42 [TARGET] = 4! 



FIG. 8: Contribution to the vertex C42 from the TARGET graph. 



Finally we calculate the contributions to the vertex N 32 - We differentiate with respect 
to U and remove the attached legs, and then differentiate with respect to V. The result is 
shown in Fig. 9. 



N 32 [TARGET] = 6 



FIG. 9: Contribution to the vertex N 32 from the TARGET graph. 



The calculation for each diagram in Fig. 1 is done in the same way as for the target 
diagram. Each graph contains a numerical factor that corresponds to the symmetry factor 
of the diagram, times the permutation factor for the left legs and right legs. When a vertex 
is substituted into one of the integral equations (21) or (22), the left legs become internal 
ones. To present the results so that the numerical factors have a more recognisable form, 
we divide each vertex by the permutation factor for the left legs. Using this notation, the 
complete results for each of the vertices are shown in Figs. 10 - 15. The results have been 
compactified by using the equation represented in Fig. 16 to remove the bare vertex Uq. 
This substitution results in the cancellation of a large number of graphs. 
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EIGHT EGG BBALL (I BBALL HAIR MERCEDES 




LOOPY LOOPY TWISTED EYE EYE EYE 



+ 




EYE EYE TARGET TARGET TARGET TARGET 



FIG. 10: The shaded box is the vertex C22 which is part of the kernel of the f2 integral 
equation. The last TARGET graph and the last two EYE diagrams should each be drawn 
as two graphs where the second one is the left-right inversion of the one in the figure. The 
graphs are combined to simplify the figure. 




FIG. 11: The shaded box is the vertex C32 which is part of the kernel of the f2 integral 
equation. The vertex C23 which appears in the \l/ integral equation is obtained by inverting 
the diagram left to right. 




FIG. 12: The shaded box is the vertex C33 which is part of the kernel of the \P integral 
equation. 
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HAIR MERCEDES TARGET EYE TWISTED 



FIG. 13: The shaded box is the vertex N22 which is part of the kernel of the \1/ integral 
equation. The function N 22 which appears in the kernel of the integral equation is given 
by the first 2 diagrams on the right side, with a full vertex on the first diagram. 



C42 



1 

U 



+ 



FIG. 14: The shaded box is the vertex C42 which is part of the kernel of the integral 
equation. The vertex C24 which appears in the integral equation is obtained by inverting 
the diagram left to right. 




TARGET EYE 



FIG. 15: The shaded box is the vertex N32 which is part of the kernel of the integral 
equation. The vertex N23 which appears in the integral equation is obtained by inverting 
the diagram left to right. 



INTEGRAL EQUATIONS FOR VARIATIONAL VERTEX FUNCTIONS 

The integral equations derived in the previous section for the vertex functions Q and \I/ 
depend on the vertices U and V. These vertices satisfy the integral equations shown in Fig. 
16, which are obtained from the equations of motion (4). In this figure, we have combined 
diagrams that correspond to permutations of external legs. 
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EGG EGG, HAIR MERCEDES EYE TARGET TWISTED 



X-XXXM'^ 

BBALL BBALLn TARGET EYE LOOPY 

FIG. 16: Integral equations for the vertices U and V. 

CONCLUSION 

In this paper we have calculated the 4-loop 4PI effective action (Eq. (3) and Fig. 1) for 
a scalar theory with cubic and quartic interactions, with a non-vanishing field expectation 
value. We derive a set of coupled integral equations that give the corresponding re-summed 
2-point vertex function. The resulting expression has the same form as the Schwinger Dyson 
equation (Eq. (18) and Fig. 3). The Kubo formulae relate transport coefficients to 2-point 
correlators. A next-to-leading order contribution to the shear viscosity could be obtained by 
solving the set of equations derived in this paper. We have checked that the equations de- 
rived in this paper are correct by expanding the re-summed 2-point vertex function to 3-loop 
order and verifying that all contributions are present, with the correct symmetry factors [38]. 

Acknowledgements: The authors thank Julien Serreau for many useful comments and 
suggestions. 



In this appendix we present some details of our derivation of the 4PI effective action for 
a scalar theory with cubic and quartic interactions, with a non-vanishing field expectation 
value. We follow the method of Ref. [18]. We use a symbolic notation that suppresses indices, 
and therefore combines terms that correspond to permutations of indices. For example, a 
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sum of three terms that corresponds to the symmetric combination of the product of a 
propagator and a field expectation value is written: 

Dij(j) k + Djkfa + Dti^j -)■ 3D (f) . (33) 

We start with the generating functional for Green's functions in the presence of quadratic, 
cubic and quartic source terms which is given by: 



Z[J, R, R 3 , R A ] = Exp [iW[J, R, i? 3 , R±] 

= j D$Exp[i(sM + J^ + ^R^ + ^R 3 ^ + ^Ra^)]. (34) 

Derivatives of the generating functional W with respect to sources produce connected 2-, 3- 
and 4-point functions which, in turn, can be expressed in terms of the proper 3- and 4-point 
vertices. 

In this appendix, we introduce some new notation for the proper vertices. First, we 
remind the reader that the definitions in Eqs. (1) and (2) were chosen to make figures 
look as simple as possible: lines, and intersections of lines, correspond directly to these 
definitions for the propagator and vertices, with no additional factors of plus or minus i. In 
this appendix we use a different definition so that the Legendre transformation has a simpler 
form. The new definition is obtained from the old one by multiplying by i, and denoted 
with a boldface character: 

{U OD , U , V , U, V} = {iU oa , iU , iV , iU, iV} . (35) 

This notation is illustrated for the bare 3-point vertex in Fig. 17. 

Uoo = i^ju^f = ^> => U 00 = i U 00 




FIG. 17: Illustration of the notation defined in (35) which will be used for all vertices in 
Appendix A. 



We will use some additional vertices that are defined in terms of the boldfaced ones in Eq. 
(35) (see Eqs (40)). These additional vertices are used only to organize the calculation in 
this appendix and never appear in the body of the article, so we do not need to boldface 
them. 
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We take functional derivatives of the generating functional W defined in (34) with respect 
to the sources, and express the result in terms of the proper vertices defined in (35). We 
obtain: 

SW 



SJ 



^4(-^ U + 3 ^ + 3), 

^ = 2-(-iD*V - 3,D 5 U 2 - iD 3 U<j) + 3D 2 + QD(f) 2 + 4 ) . (36) 
0R4 4! 

The 4PI effective action is the Legendre transform of the generating function W[J, R, R3, R4} 

5W 5W 5W 5W 
T[<j>, D, U, V] = W[J, R, R 3 , R 4 ] --JjJ-Jn R - Jr\ R z ~ Jr\ R ± ■ ( 37 ) 

It is not necessary to perform the Legendre transforms in (37) all at once: it is possible 
to implement them successively, expressing higher-order effective actions in terms of lower 
order ones. The procedure is as follows. We start with the 2PI effective action for a scalar 
theory with cubic and quartic interactions, with a non-vanishing field expectation value: 

r[0, D] = S[<j>] + ^TrLnZT 1 + ^Tr [(L>? 2 (0))- 1 {D 21 - D° 21 (<f>))} + T 2 [<f>,D] . (38) 

Since we are working to 4-loop order, T 2 [4>, D] contains all 2PI diagrams with 2, 3 or 4- 
loops. We note that the 0-dependence of ^[0,-0] can be written as a function of the 
effective classical vertex defined in Section . In the notation used in this appendix we have: 
U := U 00 + 0V O . Diagrams that contribute to T 2 [<j), D] have vertices given by U and V , 
and lines are the self consistent propagator D. 

The strategy of the calculation is to define a 2PI effective action with a modified inter- 
action: 

5W 5W 

D] = W[J, R, R 3 , R 4 ] - —J - —R . (39) 

The subscripts on Tjjy [<p, D] indicate that it depends on the sources R3 and R4 only through 
the modified interaction vertices which are defined as: 

U := U co -R 3 , V:=V -Ra- (40) 

The structure of (39) is a consequence of the fact that the cubic and quartic source terms in 
Eq. (34) can be combined with the 3- and 4-point interaction terms in Eq. (1) by making 
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the definitions in Eq. (40). The effective action defined in (39) has exactly the same form as 
the effective action in (38), with the vertices U QO and Vo replaced by U and V, respectively. 

The next step is to express the 4PI effective action in (37) completely in terms of the 
vertices U 00 , Vo, U, V and the 2PI effective action T(j V [(j), D]. This is accomplished using 
the relations: 



5T 0v _ 6W_ 5T 0v _ 5W 
5R 3 ~JR~ 3 , 5R 4 ~ 5R~ 4 



which are obtained from (39). Using the definitions (40) we obtain: 



5W 5W 
r 0, D, U, V] = T d y[^ D] - — - R 3 - — - R, 

oK 3 0/14 



ST f 



D] , T -r 



sv 



(41) 



(42) 



The last step is to to express the modified vertices U and V in terms of U and V. This 
is done by comparing two expressions that relate these vertices. 

On one hand, using (36) and (41) one obtains the expressions for the derivatives of the 
modified 2PI effective action T^y[(j),D]: 



5T 



D] 



5U 

sr D y[<f>,D) 



6V 



._(_, jD 3 u + 3jD0 + 3 ); 

— (-iD 4 V - 3,D 5 U 2 - iD 3 U<p + 3D 2 + QDcj) 2 + 4 ) 



(43) 



On the other hand, one can explicitly take derivatives of the 4-loop 2PI effective action 
with bare vertices replaced by the modified interaction vertices (Eq. (40)). The non-4PI 
contributions to the interacting part are drawn in Fig. 18 (see Ref. [39]). 



r 2 = r 4 + £ 







FIG. 18: The interacting part of the modified 2PI effective action Yfjy[(f), D] to 4-loop order. 
The 3-point vertices are U := U + <pV and the 4-point vertices are V. The 4PI terms are 
denoted T^. They are shown in Fig. 1 using the vertices defined in (2) and are not redrawn 
here. 
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Taking derivatives of the 4-loop expansion of Tjjy we obtain: 

ff w [0, g]= _g_i i i 

SU 6 2^6 4 

8 6 2 12 

gWfcgL _1 4 _ _ \ D 2 ± D ,y 1^2 1^3^ 

24^ 4 r 8 24 16 6 r 

+ - W 2 !7 + -D 5 u 2 - W 2 - -^W 3 

4 8 8 8 6 

--D 8 V(fU 3 - -D 8 U* + -D 9 <fU 5 . (44) 
2 8 12 

Equating (43) and (44) and using an iterative procedure to rearrange, we get the results in 
Eq. (45), where the dots indicate that we have dropped contributions that are of higher 
loop order. 

3? 

V = V + -D 2 V 2 + 6D 3 VU 2 - 3iD 4 U 4 + ■■■ 

3? 3 
U = U + — D 2 UV + £> 3 U 3 - -D 4 UV 2 + 6iD 5 V 3 V + 4L> 6 U 5 + • • • (45) 

Substituting (43) and (45) into (42) and using (38) we obtain from a straightforward calcu- 
lation: 

r[<f>,D,U,V\= S[<j>] + ^TrLn/T 1 + '-Tr [(D^T 1 {D*i ~ D° 2 M)} 

+r°[0, D, U, V] + T int [D, U, V] , 
r°[0, D, U, V]= ~D 2 V + ^ 3 U U + ^^ 4 VV + i,D 5 VoU 2 , 

F int [D, U, V]= -— D 3 U 2 - — D 4 V 2 - —D 6 U 4 + — D e U 3 - -D 7 V 2 U 2 
L J 12 48 24 48 8 

_I L )8 U 4 V+ ^ L) 9 U 6 (46) 

8 72 

Converting the notation for the vertices using Eq. (35), these results for r°[</>, D, V 3 , V4] + 
r m *[D, V 3 , V4] are shown in Fig. 1. Thus we find that the Legendre transform has removed 
the non-4PI terms that were present in the modified 2PI effective action (39). 
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